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Abstract 

By use of the threshold expansion we develop an algorithm for analytical evalu- 
ation, within dimensional regularization, of arbitrary terms in the expansion of the 
(two- loop) sunset diagram with general masses mj, m,2 and 7713 near its threshold, 
i.e. in any given order in the difference between the external momentum squared 
and its threshold value, (mi + m,2 + m^) 2 . In particular, this algorithm includes 
an explicit recurrence procedure to analytically calculate sunset diagrams with 
arbitrary integer powers of propagators at the threshold. 
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1. Introduction 



The purpose of our paper is to apply explicit prescriptions [|| recently obtained for 
the expansion near threshold, i.e. in powers of the difference between the external mo- 
mentum squared and its threshold value, to calculation of the sunset diagram (see Fig. 1) 
with general masses and powers of propagators. We thereby extend, to the threshold 
expansion, analysis of two-loop Feynman integrals performed in a series of papers |2|-[5] 
and based on explicit formulae for the asymptotic expansion of Feynman diagrams in 
various off-shell limits of momenta and masses |§ (see also earlier papers and operator 
analogues in ||; for brief reviews, see in H), as well as for some typically Minkowskian 
on-shell limits [ iPf . In particular, in ref. || the threshold behaviour at small (as com- 
pared to some masses not involved in the cut) non-zero thresholds was examined within 
the large mass expansion, and the description of three-particle thresholds depended cru- 
cially on a possibility to describe the threshold behaviour of the sunset diagram and a 
similar diagram with four propagators. 

The sunset diagram (Fig. 1) represents the simplest example of a diagram involving a 
three-particle threshold, at k 2 = (mi +m 2 +m 3 ) 2 . When one or two internal particles are 



massless, the (four-dimensional) results can be obtained in terms of dilogarithms [|XT| , [12|] | 
The situation gets more complicated when all three virtual particles involved in the cut 
are massive. Such a situation occurs e.g. in the two-loop off-shell contributions to the 
Higgs self energy. Although in three dimensions the result for this diagram (even with 
different masses) is quite simple |14| (cf. also in ref. |H|]), in four dimensions no exact 
results in terms of known functions (like polylogarithms, etc.) are available in this 
(totally massive) case. Moreover, there are arguments |16| that the result cannot be 
expressed in terms of polylogarithms, except for the special cases like the threshold (and 
pseudothresholds) . 

One possibility is to use various integral representations jl^, |l7Hll| , in order to get 
numerical values for given masses and the momentum. Another possibility is to study 
analytic expansions in different regions. For instance, using the algorithms presented 
in one can construct several terms of the small momentum expansion and the large 
momentum expansion. Moreover, in ref. [2(| a closed formula for the corresponding coef- 
ficients (in the case of the sunset diagram) was derived, and the occurring hypergeometric 
series were identified as Lauricella functions of three variables. Nevertheless, because of 
difficulties in constructing analytic continuation of these functions, one cannot extract 



much information about the threshold behaviour. We also note that in ref. |2T[ a differ- 
ential equation for the sunset diagram was constructed, whose derivation was essentially 
based on the integration by parts [22| and a recurrence procedure described in | 23| . The 



first step towards explicit evaluation of the coefficients of the threshold expansion was 
done in ref. ||24j| , where the threshold (and pseudothreshold) values of the sunset diagram 
(with unit powers of the propagators, and also in the case when one of the powers is 
equal to two) were calculated analytically in terms of dilogarithms of the mass ratios. 
The sunset diagram with general powers of propagators and masses is given by the 

3 Some other results for two-loop two-point diagrams with masses can be found in refs. [Of. 
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following two-loop integral (cf. Fig. 1): 

j, \= [ [ d>! d> 2 

1 ' x ' 2 ' 3j - y y [p? - mir [pi - mip k* - ft - P2 ) 2 - mir ' 1 J 

where n = 4 — 2e is the space-time dimension, in the framework of dimensional regular- 



ization [251. We shall use the notation 



(j, T = m 1 + m 2 + m 3 , f * = — , (2) 

with J2i=i 6 — 1; an d introduce the expansion parameter 

y = H 2 T -k 2 . (3) 

To make the dependence on y manifest, it is convenient to shift the loop momenta in 
eq. (HD as p\ — > pi + p 2 — ► P2 + £2^- Thus the r.h.s. of eq. ([]]) becomes 

d>i d> 2 

[p? + 2£ 1 (*p 1 ) - ei 2 2/]^ [p| + 26(^2) - HhP l(Pi + P2) 2 - 2Uk, Pi + Pa) - el2/]" 3 ' 

(4) 

Without loss of generality, in the time-like region we can choose a frame with the 
external momentum k = (k , 0) with k > 0. According to the general prescriptions of 
the threshold expansion [|l|], one should consider every loop momentum p^ (i = 1,2) to 
be of the following four types: 

hard (h): p i0 ~ \k\, \p { \ ~ \k\, 

soft (s): p i0 ~ y/y, \p % \ ~ y/y, 

potential (p): p i0 ~ y/|fc|, |pi| ~ ^/y, 

ultrasoft (us): p i0 ~ j//|A;|, |p*i| ~ y/|fc| , 

where \k\ = Vk 2 = k . 

We shall consider the most complicated case when all three masses are non-zero. 
Then only (h-h) and (p-p) regions contribute to the threshold expansion because, for 
any other region, one obtains scaleless integrals which are naturally put to zero (they 
are analogous to massless tadpoles in dimensional regularization |25|). Note that if one 
or two masses were zero then we would have non-zero (us-p) or (us-us) contributions, 
instead of the (p-p) one (cf. in refs. []I|[2"6|1). 

As the main example, we shall treat the diagram with v\ — 1/2 = 1^3 = 1. Within 
threshold expansion, this "master" sunset diagram is represented as 

L(4 - Is- 1, 1, 1) = -tt 4 - 2£ r(l + Is) k 2 £ (fey \cf- h \k 2 )- 2 ^ + Cj p - V 2e l , (5) 



j=0 



where we have extracted the factor m n ^ 2 per loop. Our results can be easily generalized 
for sunset diagrams with any indices fj. 
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The (h-h) contribution is given by expanding the integrand in (B) in Taylor series in 
y, so that these are sunset integrals, with various (higher) indices of the lines, evaluated 
exactly at the threshold (y = 0). In sections 2-4 we shall describe how an arbitrary 
integral of this kind can be analytically evaluated. It happens that general solutions of 
recurrence relations constructed in |23[] cannot be directly applied at the threshold. Nev- 
ertheless, those relations are useful at the threshold, providing equations which should 
be used together with the usual integration- by-parts relations [22|. 



The (p-p) contribution is given by expanding the propagators in (|J) in Taylor series 
in the squares of time components, pf > P20 an d (pio + P2o) 2 - The evaluation of the (p-p) 
part of the expansion is presented in Section 5. In Section 6 we compare our results 
for the threshold expansion with results based on numerical integration. We briefly 
summarize our results in Section 7. 

In general, the threshold expansion has a complicated structure and it is difficult to 
see how it works. Nevertheless, in our particular example we can present some simple ar- 
guments. Namely, it happens that the remainder of the expansion of the sunset diagram 
can be constructed by applying, to eq. (f|), an operator 

Here Ty- N ' denotes an operator that picks up the first N terms of Taylor expansion in y 
(cf. eq. (H) below), whereas corresponds to the first N terms of the expansion in 

PiQ 

squares of time components pf Q , p\\ and (P10+P20) 2 (cf- ec L- (E3) below). It is assumed that 
these operators "commute" with the the loop integrations, i.e. they should be applied 
directly to the denominators in eq. (f§). 

For instance, for the integral with ^ = 1 we can write 

L(n; 1, 1, 1) = (l - 7Z {N) ) L(n; 1, 1, 1) + TZ (N) L(n; 1, 1, 1) 

= ( T v N) + ~ T v iN)T t 2) ) ^ + X > V ' W 

where the terms Tj-^L and T^~ 2 ^L give us the (h-h) and (p-p) contributions to the 
threshold expansion (|5|), respectively. For dimensionally-regularized integrals, the prod- 
uct of these operators, vanishes. To see this, one can integrate in p w and 

y PiO 

P20 using residual theorem: then, the resulting integrals in pi and P2 are scaleless (i.e. 
analogous to massless tadpoles, but in (n — 1) dimensions), and therefore vanish. 

Finally, the last term on the r.h.s. of eq. (0), IZ^Lin; 1, 1, 1), plays the role of the 
remainder of the expansion. To show this, consider first the region of small loop momenta 
where the operator T^ N ^ is dangerous because it produces infrared threshold singularities. 

However these singularities are removed by the operator f 1 — 2 ^ which increases 
powers of the variables p\ and P2 in the numerator. Then, consider the region of large 
loop momenta where the operator ^ is dangerous because it generates specific 
divergences (corresponding to ultraviolet singularities in lower dimensions). Still, the 
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latter are removed by the operator (l — 1y- N n which effectively increases the powers of 

integration momenta in the denominators. Thus the remainder TZ^ does not involve 
new divergences, neither ultraviolet nor infrared ones, as compared to the initial Feynman 
integral Moreover, taking into account the operator (l — Ty N ^ we see that the order 
of this remainder is y +1 , up to logarithms. 

2. Taylor expansion and the (h-h) contribution 

The formal expansion in y of the denominators in the integrand of eq. @ can be 
performed via 

1 x (//,). rf i- yh 

[pf + 2&fe) - &vr = £ if \pi + 2Uk P i)r +M ' (8) 

where = r(u + j)/T{y) is the Pochhammer symbol. If we denote P3 = —p± —p2, the 
third denominator can also be taken into account by eq. (§). Then, collecting the terms 
with given powers of y, we get the (h-h) contributions (cf. eq. (|5|)). According to the 
r.h.s. of eq. each separate term corresponds to the integral (f|) with shifted powers 
of propagators (z/j — > Vi and y = 0, i.e. at the threshold. 

One should however keep in mind which variables are considered as independent ones. 
Let us consider £j, eq. (0), as "external" parameters characterizing the mass ratios. Then 
the remaining variables are y, k 2 and /i^ = (mi + ni2 + rriz) 2 , which are dependent, 
according to the definition (H) . Since y is the main parameter of the threshold expansion 
(cf. eq. (^|)), we can choose between the following pairs of independent variables: (i) y 
and k 2 , and (ii) y and /z^. 

On one hand, in the context of the threshold expansion (including the related 
issue of the non-relativistic limit), the first set, y and k 2 , looks more natural (cf. eq. (0)). 
In particular, the momentum integrals arising after expansion in y (H) are functions of 
k 2 . In section 5, we shall also see that such a choice is more convenient for evaluating 
the (p-p) contribution. On the other hand, when evaluating the threshold values of the 
momentum-space integrals (corresponding to the (h-h) contribution), it is convenient to 
put k 2 = i.e. to calculate 

L T (n; Vl , z/ 2 , u 3 ) = L(n; v x , »2, ^)| fe 2 =(mi+m2+m3) 2 . (9) 

This would correspond to the second set of independent variables, y and /i^. In partic- 
ular, when v; = 1 we define 



y - J 



L(4 - 2e- 1,1,1) = -^T 2 {1 + e) & ^ (J-J ^ h ~ h) (/4) + ^"V 



• (10) 

In such way, it is more straightforward to use the threshold results for the lowest integrals 



obtained in pl| . Another advantage, as we shall see, is connected with the "book 



keeping" of ultraviolet singularities. 
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Anyway, the transition between these two options (and, in particular, between the 
coefficients Cj and Cj in eqs. (]5|) and (|10D, respectively) is rather straightforward. For 
instance, any integral Lt (0) can be presented as (/if-)™ ^ " 2 ^ times a function of 
dimensionless parameters £j, eq. @. When we consider the momentum integrals corre- 
sponding to the expansion (j|), the only difference is that we need to replace ^ by k 2 , 
namely: 



d>i d> 2 



\p\ + 2UkpiT [pi + 2!; 2 (k P2 )r [( Pl +p 2 ) 2 - 2Uk,Pi +P2T 

= (fc 2 //4)"-^-^ L T {n-u u u 2 ,u 3 ) . 



(11) 



If we consider the (h-h) contributions (§) corresponding to L(n; 1, 1, 1), the first terms 
are proportional to 

L T (n; 1,1,1), 

m\L T (n] 2, 1, 1) + m\L T (n] l,2,l)+mlL(n;l,l,2), 

m^m^LT^] 2, 2, 1) + m^mgLT^n; 1, 2, 2) + m\m\LT{n\ 2, 1, 2) 

+r<L T (n; 3, 1, 1) + m\L T (n; 1, 3, 1) + m\h T {n; 1, 1, 3), (12) 

etc. It is easy to see that each new order will contain ultraviolet-divergent combinations 

ml j L T (n;j + 1, 1, 1) + m 2 2 j L T (n; l,j + 1, 1) + m^L r (n; 1, 1, j + 1). (13) 

The situation is different if we use the variables y and /x^, according to eq. ([T0|). 
It is instructive to consider the formal Taylor expansion around the threshold k 2 = n\ 
(» = 0), 



i=o 



2/V 
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<% 2 



(14) 



L{n\v u v 2 ,v 3 ) 

k 2 ={mi+m2+m^) 2 

which also corresponds to the (h-h) contribution. An algorithmically convenient way to 
construct expressions for the derivatives in ([14]) is to use the formulae 



'A 

dk 2 



L(n; z/i, u 2 , v 3 ) 



■IV 



7T 



2.7 



>i)i(^)i(^)i ^(n + 2j; z/i + j, v 2 + j, 1/3 + j) (15) 



and 

7r 2 A; 2 (i/i + j — 1)(^ 2 + J - l)(^s + J - l)^(n + 2j; ^ + j, v 2 + j, i/ 3 + j) 

= (vi + 1/ 2 + ^3 - n + j - 1) L(n + 2j - 2; v x + J - 1, v 2 + j - 1, u 3 + j - 1) 
+ (z/i + j - l)m 2 L(n + 2j - 2; ^ + j, u 2 +j- 1, 1/3 + 3 ~ 1) 
+ (^2 + j - l)ml L(n + 2j -2;vi+j- 1, z/ 2 + j, 1/3 + J - 1) 
+ (ia, + j — l)m 2 L(n + 2j-2;v 1 +j-l,v 2 +j- 1, z/ 3 + j). (16) 
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Applying eq. ([16]) j times, we reduce the space-time dimension back to n = 4 — 2e. Both 
formulae (|T5|) and flUf ) can be easily derived by using the modified Feynman parametric 
representation for the integral (|l]) given by eq. (4) of ref. |23| (cf. also in refs. [27]). 
When we consider the Taylor expansion of L(n; 1, 1, 1), the derivatives are given by 

-k 2 -^-L{n- 1, 1, 1) = (3-™)L(n; 1, 1, l)+m 2 L(n; 2, 1, l)+m 2 L(n; 1, 2, l)+m 2 L(n; 1, 1, 2), 
ok z 

(17) 



9 \ 
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\(k 2 ) 2 L(n; 1, 1, 1) = |(4 - n)(3 - n)L(n; 1, 1, 1) 

+(4-n) [miL(n; 2, 1, 1)+ m 2 2 L(n; 1, 2, 1) + m\L(n\ 1, 1, 2 

+m 2 m 2 L(n; 2, 2, 1)+ m 2 m 2 L(n; 1, 2, 2) + m 2 m 2 L(n; 2, 1, 2) 
+m^L(n; 3, 1, 1) + m\L(n- 1, 3, 1) + m\L{n; 1, 1, 3), (18) 

etc. In these derivatives (see eq. (|14T)), we should put k 2 = (mi + m 2 + m 3 ) 2 (L — > L T ). 

Note that now the ultraviolet singularities are present only in the y° and y 1 terms 
of the expansion of L(n; 1, 1, 1). One can see that, apart from the lower integrals, the 
derivatives flTTD and (|1^) contain the same combinations as (|T2|). Furthermore, one can 
get the results (|17D and (|T8|) directly from eq. (p!2|), by multiplying each contribution by 
the factor (|TT|), (k 2 //j^) n ~ 3 = (1 — y/ ^j) n 3 , expanding it in and collecting terms at 
given powers of y. 

The analytic results for L r (4 — 2e; 1, 1, 1) and L T (4 — 2e; 1, 1, 2), expanded in e up 
to the finite part, are given in eqs. (23) and (22) of ref. [p^j , respectively. The results 
for L T (A - 2e; 2, 1, 1) and L T (4 - 2e; 1, 2, 1) can be obtained from L T (4 - 2e; 1, 1, 2) via 
permutation of the masses ??v Those results involve the following functions: 

T~ (z) = Li 2 (-z) - Li 2 (-l/z) + 2 In z ln(l + - In 2 z 
= 2Li 2 {-z) + |vr 2 + 2 In z ln(l + z) - \ In 2 z 
= 2Li 2 (1/(1 + z)) - |vr 2 + ln 2 (l + z) - \ In 2 z, (19) 



^ = arctan(mj mi + m2 + m3 l, 9 1 + 9 2 + 9 3 = n. (20) 
\ V mim 2 m 3 / 

Note that the arguments of T~ are just the mass ratios mj/mi, and the following inversion 
property is valid: T~(mj/mi) — —T~ (mi /rrij) . In particular, T - (l) = (remember that 
Li 2 (—1) = — j7>7r 2 ). Similar functions have also appeared in refs. p8| . 

In the next two sections, we shall consider evaluation of the integrals Lt with higher 
powers of propagators. 
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3. Recurrence relations at the threshold 



In this section, we shall consider relations connecting the integrals (|l|) with differ- 
ent indices z/j. When L is used without its arguments, the non-shifted n and Ui are 
understood. The standard notation for the raising and lowering operators reads 

1+L = L(n; Vi±l, u 2 , u 3 ), 2+L = L(n; u u z/ 2 ±l, u 3 ), 3+L = L(n; Ui, u 2 , v 3 ±l). (21) 



When the external momentum squared takes its threshold value, k 2 
we shall use the notation ©. 
Using the identities |22| 



// 







A (J) 

A 4 



dpia [ [p? - ml}" 1 [p| - m lY 2 l(k -p%- p 2 ) 2 - m 2 ] 



2] u 3 



(mi + m 2 + m 3 ; 



0, (22) 



with A^} = \px fM ,P2 l j,, k^, we obtain a set of integration-by-parts relations. For example, 
one of the equations can be presented as 



2m 2 1 u 1 l + + (ml + m 2 2 + ml- k 2 )v 2 2 + + 2m^ 3 3 + + 2m ^ lV2 3-1+2 

v% — 1 



2n 



v 3 

, n — 2 1/1 — v-i + 1 

2z/x - i/ 2 - 2z/ 3 - z/ 2 l"2 + + v 2 2TT 

v 3 -l 



L 
L. 



(23) 



and five other equations can be obtained from (p3|) by permutations of the indices (1,2,3). 

Using three of these six equations, we can e.g. express 3~l + 2 + L, 2 _ l + 3 + L and 
l~2 + 3 + L in terms of 1 + L, 2 + L and 3 + L. However, three remaining equations for 1 + L, 
2 + L and 3 + L happen to be linearly dependent, for arbitrary k 2 . Useful corollaries of 
eq. (p3|) and its permutations are 



m?z/ 1 (z/ 1 + l)l + l+ - m 2 3 v 3 (is 3 + 1)3+3+ 
m 2 i z/ 2 (z/ 2 + l)2+2+ -m^ 3 (zy 3 + 1)3+3+ 



L = \ (n-2vi-2)z/il+ - (n-2v 3 -2)v 3 3 + L, 

(24) 

L = \ \(n-2v 2 -2)v 2 2 + - (n-2u 3 -2)u 3 3 + ] L. 

(25) 

Below we shall use these two equations as a part of the main algorithm. 

Let us denote a = v\ + v 2 + v 3 . Since the threshold values of the integrals with a = 3 
(vi = v 2 = v 3 = 1) and a — 4 [yx = v 2 = 1, v 3 — 2 and permutations) are known ||24j| , the 
first non-trivial step is to get results for six integrals with o = 5 (L(n; 2, 2, 1), L(n; 1, 1, 3) 
and permutations). Introduce the following notation for the symmetric sums of these 
integrals: 

S 221 = L{n; 1, 2, 2) + L{n; 2, 1, 2) + L{n; 2, 2, 1), (26) 
S U 3 = m\L(n\ 3, 1, 1) + m 2 2 L(n; 1, 3, 1) + m 2 3 L(n; 1,1,3). (27) 
Then, using eqs. (|2~4])-([2"5]) we get 



L(n; 1, 1, 3) = -L {5 113 + 1 (n - 4) [2L(n; 1, 1, 2) - Lin- 2, 1, 1) - L(n; 1, 2, 1)]} (28) 
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and analogous results for L(n; 3, 1, 1) and L(n; 1, 3, 1). Furthermore, using eq. (p3|) and 
its permutations, we get 



L(n; 2, 2,1) 



2 

mf 



m 2 + m 3 



{277^221 + |^ii3 - \{n - 4)L(n; 1, 1, 2) 



-\{An - 13) [L(n; 2, 1, 1) + L(n; 1, 2, 1)] + L(n; 2, 2, 0) J 



(29) 



and analogous results for L(n; 1, 2, 2) and L(n; 2, 1, 2). Note that at the threshold 



(k 2 — m\ — m\ + mjj 



fc 2 = (mi +m2 ) 2 



2(m x + m 3 )(m 2 + m 3 ). 



Furthermore, taking the sum of eq. (^) and its permutations, we can express S 22 \ 
via an d lower integrals. In particular, at the threshold we get 



113 



^(5n-17)/x T [L T (n; 2, 1, l) + L T (n; 1, 2, l)+L T (n; 1, 1, 2)] 



-\{n - 3) [miL T (n; 2, 1, 1) + m 2 L T (n; 1, 2, 1) + m 3 L T (n; 1, 1, 2)] 
+| [(ma + m 3 )L(n; 0, 2, 2) + (mi + m 3 )L(n; 2, 0, 2) + (m x + m 2 )L(n; 2, 2, 0)] , (30) 

where the superscript "(T)" means that the corresponding sums, S221 and ^ji^, are 
considered at the threshold. Therefore, using the integration-by-parts relations we reduce 
all the six integrals with a = 5 to a single unknown function (say, S221 )• 

A missing "block" completing the algorithm can be obtained using some equations 
presented in ref. [ 23| . Consider eqs. (74) and (81) of p3| . Their l.h.s.'s contain a factor 



-D123 = 



k 2 — (mi + ^2 + m 3) 2 k 2 — (—mi + m 2 + m z) 2 
k 2 — (mi — m 2 + m 3 ) 2 k 2 — (nil + m 2 — m 3 ) 2 



(31) 



In the case of interest (i.e. at the threshold) these l.h.s.'s vanish, since -D123 = 0. There- 
fore, at the threshold these equations cannot be used in the way as it is suggested in 
ref. |23j, since this is a degenerate case. Nevertheless, the fact that the r.h.s.'s should 
also be equal to zero, provides some non-trivial conditions on the integrals involved. 
It is interesting that at the threshold both eqs. (74) and (81) of [23] lead to the same 
condition: 

\h>inii \(2n — u 1 — 2u 2 - 2z/ 3 - l)mi + (n - v 2 - 2u 3 )m 2 + (n - 2u 2 - z/ 3 )m 3 ] 1 + 
+v 2 m 2 [(n — ui — 2v 3 )nii + (2n - 2v x - v 2 - 2z/ 3 - l)m 2 + (n - 2v x - ^ 3 )m 3 ] 2 + 
+z/ 3 m 3 \(n — v x — 2v 2 )nii + (n — 2v x - v 2 )m 2 + (2n - 2v x - 2v 2 - u 3 - l)m 3 ] 3 + 
— z/2f3m2m 3 l~2 + 3 + — i/ii/ 3 mim 3 2 _ l + 3 + — ^iz/2mim23~l + 2 + 

-\(n -vi-v 2 - u 3 )(3n - 2v x - 2v 2 - 2v 3 - 2)} L T (n; u h u 2 , u 3 ) = 0. (32) 
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For example, when v\ = v 2 = 1/3 = 1, eq. ([32]) yields 

(n - 3) {m 1 (2m 1 +m 2 + m 3 )L T (n; 2, 1, 1) + m 2 (m 1 + 2m 2 + m 3 )L T (n; 1, 2, 1) 

+m 3 (m 1 + m 2 + 2m 3 )L T (n; 1, 1, 2) - ±(3n - 8)L T (n; 1, 1, 1)} 
= m 2 m 3 L(n; 0, 2, 2) + mim 3 L(n; 2, 0, 2) + mim 2 L(n; 2, 2, 0), (33) 
where on the r.h.s. we have just the products of massive tadpoles, 

L(n; v u u„ 0) = (m 2 )"/ 2 "- (^)»^ ^ = l5^7^ = — (34) 

and permutations. We have checked that the results for Lt(ti; 1, 1, 1), Lt(ti\ 1, 1, 2) (and 
permutations) presented in eqs. (22) and (21) of |24| satisfy eq. fl33p. 



Furthermore, considering eq. ([52]) for z/l = i/ 2 = 1, v 3 = 2 (and permutations), we 
obtain an extra (independent of (|30|)) condition on the sums S221 and S^ 3 \ In this way 
we can obtain results for the integrals L T (n; 2, 2, 1), L T (n; 1, 1, 3) (and permutations) in 
terms of lower integrals. However, at this step (solving the system of linear equations) 
we get the factor (n — 4) in the denominator. This means that we need to know the e part 
of the lower integrals, in order to get the finite part of Lx{n\ 2, 2, 1), etc. Alternatively, 
we can calculate the finite part of L^(n; 2, 2,1) via straightforward calculation. As we 
shall see, this will be enough to calculate the higher integrals (see in Appendix A for 
details). 

4. Calculation of L^(4 — 2e; 2, 2, 1) with different masses 

We start from the two-fold integral representation (5) of ref. [23], which in the case 
v\ = V2 = 2, v 3 = 1 reads 



L T (4 - 2e; 2, 2, 1) = 7r 4 ~ 2£ T(l + 2e) (mim 2 )^ mf £ 

x 11 ^ dV g ^ (35) 

^ (mi^ + m 2 r] + m 3 )^ 3£ [m 2 m 3 £(l-??) 2 + miin^l - £) 2 + m 1 m 2 ((-f;) 2 ] 

The threshold singularity originates from the region £ ~ 77 ~ 1. Indeed, analyzing eq. (5) 
of p4| we see that the threshold singularity takes place (in the four- dimensional space) 
when V\ + z/ 2 + u 3 > 5. It is of an infrared origin and (as we shall see) appears in 
dimensionally- regularized integrals as a pole in e = (4 — n)/2. 

The double integral in the second line of eq. (|35D can be decomposed as 

1 77 dzd V 



3e { I [m 2 m^(l - 7]) 2 + m 1 m 3 ?7(l - £) 2 + m 1 m 2 (^ - i]) 2 ] 1+2£ 

OO OO 

1 r r d£ dr] [m 1 (l— +m 2 (l— 77)] 

/ij- 7 7 [m 2 m 3 ^(l— r^) 2 + m 1 m 3 r](l — ^) 2 + m 1 m 2 (^ — ?]) 2 ] (m^ + m 2 r] + m 3 ) 

+0(e), (36) 

10 



where the threshold singularity is only in the first term (which is simpler than the original 
integral fl35D ), whereas the second term is finite as n — > 4. 

To evaluate the integrals in (|36|), we can proceed in the same way as in ||24)| . However, 
the following substitution of variables^ happens to be more efficient than eq. (13) in ||24j| : 

£ = 1 + (m 2 /m 3 )Az, 77 = 1 + (mi/m 3 )A(l - 2) . (37) 
Then, the region of integration in variables z, A is defined by 

— m 3 /(m 2 A) < 2 < 1 + m 3 /(mxA) (for A > 0) , 
1 + m 3 /(miA) < z < -m 3 /(m 2 A) (for A min < A < 0) , 

with A m j n = — m 3 (mi + rn.2 )/ (777.1 777.2). Note that the resulting integrand is invariant under 
(z, 777.1) (1 — 2,777,2). The threshold singularity corresponds to the region A ~ 0: the 
transformed integrand of the first integral in (|36|) contains A~ 1_4e . 

Evaluating the integrals in eq. (EB) we arrive at the following result: 



L T (4-2e;2,2,l 
x < -777,3 



7T 4 - 2e T 2 (l + e) (m 1 m 2 m 3 )- 3£ /^ 2+5£ 



mi + m 2 + m 3 



/mi 
\m 3 



mim 2 m 3 
Vm 3 



7T 



+ 



2£ 

4tt 2 



+ 7r£(mi, m 2 , "73) + C(6»i, 6> 2 , 9 3 

ah 1 2 m l 1 1 2 m l 1 1 2 m 2 

-An9 3 + hr ±m ± kf — 



m 2 



m 3 



m 3 



+ 0(e), (38) 



with 

c{e u e 2) e 3 ) = ci 2 (2e 1 )+c\ 2 (29 2 ) 

and 

£(7771,7772,7773) = In 



Cl 2 (29 3 ) + C1 2 (tt-20i) + C1 2 (tt-2# 2 ) + C1 2 (tt-2# 3 ) 

(39) 



(7771 + m 2 )(m 2 + m 3 )(m 3 + 



6 In 2. 



■ , (40) 

y (mi + m 2 + m 3 J d J 

The functions T~ and 9i are defined in eqs. ([19]) and (pD|) , respectively (cf. eqs. (14) and 
(20) of [0). The Clausen function is defined as Cl 2 (9) = Im Li 2 (e l6 j (for details, 



see e.g. in [^J). Note that the combination of Clausen functions ( p9|) is related to the 
volumes of asymptotic tetrahedra in hyperbolic space of constant curvature (see in fl30|). 
Worth noting is that at the threshold the result for the symmetric combination (EBT) is 



q(T) 
°221 



7T 



4-2e 



r 2 n 



(mim 2 m 3 ) 1 / 2+3e / u^/ 2 5e 



7T 



— + 7r£(mi, m 2 , m 3 ) + C(6»i, 6> 2 , 3 ) 



(41) 



Once more, we would like to emphasize that the 1/e poles in eqs. (^) and ([41]) (accom- 
panied by an extra factor of tt) correspond to the threshold singularity: there are no 
ultraviolet divergences in these results. 

4 A. D. is grateful to J.B. Tausk for discussion of this substitution in connection with ref. E3). 
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Using recurrence relations, we can obtain results for L^(4 — 2e; 1, 1,3) with differ- 
ent masses, as well as for the higher integrals. Some relevant results are collected in 
Appendix B. 

In the equal-mass case, 61 = 62 = 63 = (cf. eq. (f2Cf )), T _ (l) = (cf. eq. (|i~9"D), 
C(7r/3,7r/3,7r/3) = 5Cl 2 (vr/3) (we take into account that Cl 2 (2tt/3) = |C1 2 (vr/3)). 
Therefore we get 

Lr(n; 2, 2, l)| mi=TO=mjsm = Lp>(n; 2, 2, 1) 

= r2(1 +e) m — (i _ ln3 - 61n2 ) + ^Cl 2 (I)} + 0(e). (42) 

It is interesting that Cl 2 (vr/3) (i.e. the maximal value of Cl 2 (6)) also appears in the 
two-loop integrals containing a threshold singularity^. 

In fact, knowing the result (03) we can obtain e-part of the integrals L^ q \n; 1, 1, 2) 
and L^ q \n; 1, 1, 1). Using recurrence relations, we get 

4 eq) (n; 1, 1, 2) = ^— 1 — {-8m 2 (n - 4)4 q) (n; 2, 2, 1) + (2n - 7)L^ (2, 2,0)}. 

(n— 3)(3tt, — 10) L J 

(43) 

It is important that the integral (fO) enters with a factor (n — 4). Expanding in £ we get 



4 eq) (4 - 2e; 1, 1, 2) = vr 4 ^ T 2 (l + e) - ^ - 1 + I 



11 

2 ■ y £ 



l7i [l + 5e-eIn3-6eln2]--^Cl 2 (^) }> + C(e 2 ). (44) 



3^ J 3V3 U 

Furthermore, using eq. ( |3~3f ) we obtain 



4 eq) (n; 1, 1, 1) = _ 3 6 } ^ _ {4(n - 3) 4 q) (n; 1, 1, 2) - L<«0 (n; 2,2,0)}. (45) 
Finally, we arrive at 

4 q) (4 - 2e; 1, 1, 1) = n 4 ~ 2s T 2 (l + e) m 2 " 4e (- JL - - + - + — e 
T v ' ' ' ; v 1 I 2e 2 4e 8 16 

+i| [-2 - 13e + 2e In 3 + 12s In 2] - ^|ci 2 (|) \ + (9(e 2 ) . (46) 



5 For some other examples illustrating occurence of the Clausen function in two-loop calculations, 
including vacuum integrals and the pseudothreshold values, see e.g. in H,||. Note that Cl 2 (tt/3) = 
ty'(§) - §tt 2 ] /(2V3) (cf. in H|i). 
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5. Evaluation of the (p-p) contribution 

According to the prescription of [|lj, we should start by expanding the propagators 
in the integrand of eq. @ in p? , 

I =k V ^ (A7) 

where, as before, we imply that p 3 = — pi — p 2 - However, it happens inconvenient to 
expand integrand in the time components of the loop momenta (according to eq. Q47D ) 
and then take pio and P20 integrals, because one arrives at rather cumbersome integrals 
in pi and f>2- We shall proceed in a slightly different way. 

Let us instead exponentiate, by use of the alpha parameters, every propagator of the 
unexpanded Feynman integral @) and perform integration in p\ and f>2- In particular, 
for the integral L(4 — 2e; 1, 1, 1) we have 



00 00 00 



;2e n-1 



daida2do;3 



1 7T / / / 7 Cxp 



-iy(Ci«i + ^2«2 + £f « 3 ; 



000 



dpi dp 2 oexp |i[(«i + a 3 )pi + 2a 3 j9i p 2 o + («2 + "3)^20]} 

x exp {i [2/c Pio(6ai - £3^3) + 2fc P2o(6"2 - £30:3)]} • (48) 

Only now we perform Taylor expansion: namely, in the exponent with pl , p\ Q and 
PioP20- Then the evaluation of the integral of an arbitrary resulting term is simple: we 
perform integrations in pio and P20 and obtain terms proportional to derivatives of two 
5-functions, <5(2/c (£iCti — £3^3)) and 5(2A; (£ 2 tt2 ~ 6 a 3))- Then, we take the integrals in 
«i and 02, and finally in a 3 . As a result we have obtained an explicit formula for the 
(p-p) contribution in arbitrary order, written through a finite eight-fold sum. 

In eq. (||) the (p-p) contribution starts from the order y 2 , so that p ' ) = C[ p ^ = 0. 
For the first two non-trivial orders we have 

r (p-p) _. ^(^i^) 172 "' (AQ , 
° 2 ~ 4,(1 -e)(l- 2e) ' (49) 

Ct P) = ^-eXi-^) [ (1 " 2£) - (5-186)666] • (50) 

In the equal-mass case, we get 



£r(p-p),(eq) _ ( \y +1 71 
3 



a j ( — I- 3 In 3 ) + bj 



0(e) . (51) 



The coefficients aj and bj (j < 12) are given in Table 1. Note that the threshold expansion 
of the imaginary part of the sunset diagram is completely characterized by the pole part 
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Table 1: The coefficients in eqs. (|5TD, fl55|) and fl56|). 



j 


cij 


bj 


CLj = a' 


bj 


b'i 

J 


°j 


2 


j_ 

12 


l 

4 


J_ 

12 


2 
9 


l 

4 


_J_ 

12 


3 


J_ 

48 


J_ 

16 


J_ 

16 


41 
144 


_3_ 

16 


1 

8 


4 


7 
768 


15 

512 


13 
256 


1279 
4608 


79 
512 


11 

64 


5 


1 
192 


9 
512 


11 

256 


2407 
9216 


17 
128 


59 
320 


6 


83 


583 


305 


72019 


5767 


1901 


24576 


49152 


8192 


294912 


49152 


10240 


7 


233 


1693 


1073 


269359 


20719 


261607 


98304 


196608 


32768 


1179648 


196608 


1433600 


8 


1381 


6889 


7623 


1345217 


100361 


2034657 




786432 


"1 f\ A O P" *T/^ 

1048576 


262144 


6291456 


1048576 


11468800 


9 


2129 


32695 


13623 


52988741 


183713 


27534999 


1572864 


6291456 


524288 


264241152 


2097152 


160563200 


10 


108257 


1701199 


781899 


3184520135 


10819013 


1697828907 


100663296 


402653184 


33554432 


16911433728 


134217728 


10276044800 


11 


352373 


5653147 


2809445 


35912515043 


119896193 


14363946379 


402653184 


1610612736 


134217728 


202937204736 


1610612736 


90429194240 


12 


4677635 


382328867 


40375137 


900245481797 


2959291009 


220597949391 


6442450944 


128849018880 


2147483648 


5411658792960 


42949672960 


1446867107840 



of the (p-p) contribution, i.e. by the coefficients aj in Table 1. We have compared our 
results for the imaginary part with analytical expressions (in terms of elliptic integrals) 



presented in |18|] (cf. also in p4| ) and found complete agreement. 

Note that for an arbitrary, L-loop "water-melon" diagram (consisting of two vertices 
and (L + l) massive lines between them) the (p-p-. . .-p) contribution can also be analyt- 
ically calculated in every order. In particular, in the leading order we have the following 
result: 

iV n + i),/2 r((£ _ 3/2)L + 1} ft| e i/ 2 -,j (fc2)L/a ^ e _ 3/2)L+i , (52) 

where, by analogy with eq. (@), £j = rrij/(J2i=i m i)- 

6. Results and numerical comparison 

Collecting parts of the algorithm described in the previous sections, we obtain terms 
of the threshold expansion of the sunset diagram (|l|). To see how efficiently the threshold 
expansion works, we have compared our results with a semi-numerical program based 



14 



on the algorithm of [|T^]q The ultraviolet-divergent terms were compared analytically, 
whereas the finite (in e) part was treated numerically, using an integral representation 



from [19|. 

Our approximations correspond to the threshold expansion given by eqs. (|]) and (|1 
We refer to the sum of terms up to order y N as to an iV-th approximation. We illustrate 
the convergence in the plots, Fig. 2 and Fig. 3, where we present the subtracted (i.e., 
without poles in e) real part and the imaginary part of (/iy) _1+2e Lr(4 — 2e; 1, 1, 1). Note 
that the imaginary part is finite as e — > 0. These quantities are shown as functions of 
fc 2 //x|. (which equals k 2 /(9m 2 ) in the equal-mass case). In all plots, the highest-order 
curve is denoted by the solid line, whereas the lower approximations are drawn by various 
dashed lines (see the plots). The results of the above-mentioned numerical program [19 
are displayed as crosses. 

As the first example for numerical comparison, we have chosen the case £i = 0.1, 
£ 2 = 0.3, £3 = 0.6, when the masses are essentially different (but still, one cannot neglect 
either of them). Here, the expansion (|T0| ) in terms of y/ ^ is used. The results are 
shown in Fig. 2. On one hand, we can see that the expansion of the real part converges 
reasonably well. On the other hand, the expansion of the imaginary part breaks down 
just above k 2 / ^ ~ 1.5. In fact, the same happens to the real part, when k 2 j — 2.2. 
This is because the expansion parameter y/fj^ is already not so small. 

If we use, for k 2 > fi 2 ^, the expansion in terms of y/k 2 the convergence is much better. 
For instance, at k 2 = 2^ the N = 7 approximation provides a four-digit precision for 
the imaginary part. In general, it appears that an optimal choice is to use the expansion 
in y/fi^ (cf. eq. fliTf)) for k 2 < and the expansion in y/k 2 (cf. eq. (|5|)) for k 2 > 
For the imaginary part, it is enough to use only the y/k 2 expansion, since it vanishes for 
k 2 < [i 2 T . 

As the second example, let us consider eq. (|I~0"D in the case of equal masses, mi = 



m 2 = m 3 = m (& 
contributions are 



1/3), remembering that C { 



(p-p) 



(p-p) 



0. The two lowest (h-h) 



a 



(h-h),(eq) 




(h-h),(eq) 



,1 £ 



1 5 St 



Qe 2 As 8 9v/3 

+ 0(e). 



+ 0(e), 



1 23 4vr 
Te + 24 ~ 



(53) 
(54) 



It is easy to check that these two terms absorb all ultraviolet singularities of the original 
sunset integral. For j > 2, the coefficients have the following form: 



^(h-h),(eq) 



7T 

7s 



oj ( - + 31n3- 61n2) + b 3 



+ 



7T 



c 



(p-p)>(eq) 



7T 
73 



a', - + 31n3 )+K 



£ 



3 



+ 0(e) 



^Ch^j -c j + 0(s), (55) 

(56) 



6 We are grateful to P. Post and J.B. Tausk for kind permission to use their REDUCE program for 
the comparison. 
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The 1/e poles in the (h-h) contribution, eq. (|55|), correspond to the threshold singu- 
larity. However, since there is no threshold singularity in the case z/ x = z/ 2 = z/ 3 = 1, these 
poles should cancel with those from the (p-p) contribution, eq. (0). This is the case if 
a'j = a,j. We have checked that this property is valid for all available contributions. The 
corresponding coefficients (for j < 12) are collected in Table 1. 

As a result, using eq. fllPf ) we get the following threshold expansion in the equal-mass 
case: 



L (eq) (4 - 2e- 1, 1, 1) = -7r 4 - 2e r 2 (l + e)(9m 



2 , J -2 : J £t(h-h),(eq) + V ^(h-h),(eq) 



2 

V3 



^ ln 9^-^ ln2 + 5Cl2 U 



.4 



f y 



\9m 2 



71 



v/3 V9m 2 



9m 2 



C 



9m 2 



with the "form factors" 



■0(e),(57) 



^(*) = E 



Eft 



c{z) = E 



(58) 



Note that the imaginary part of the expression in braces in fl5"T|) equals 

2tt 2 



v/3 V9m 2 



^(fc 2 - 9m 2 



(59) 



where the function A(z) can be expressed in terms of elliptic integrals JT8 



We have also compared this equal-mass example with the program based on the 



algorithm of [p^| . The results of this numerical comparison are shown in Fig. 3. We have 
used the "combined" expansion: below the threshold (for k 2 < /i 2 ^) our approximations 
correspond to the expansion (|5TD-(|58|), whereas beyond the threshold (for k 2 > fxj) we 
switch to the expansion in y/k 2 . Technically, this transformation can be done just by 
substituting the arguments of the functions (|58|) via y/{jj^) — > (y/k 2 )/(l + y/k 2 ) and 
re-expanding in y/k 2 , up to the given order N. 

From the plots it is clear that subsequent approximations are getting better in a 
wide range of the values of the external momentum squared, i.e. the threshold expansion 
indeed works well. Close to the threshold, it is enough to take just a few terms of the 
expansion for achieving precision which is better than that of the numerical program . 
For k 2 = 0.5/i 2 -, the N = 7 and iV = 12 approximations to the real part reproduce four 
and six digits, respectively. When k 2 = 2/i 2 ^, the N = 6 approximation gives us three 
digits, whereas the iV = 12 curve provides six-digit precision. 

Note that the threshold expansion (at least, in its present form) does not seem to 
work in Euclidean region, i.e. when k 2 < 0. Then, the limits when some of the particles 
become massless are also non-regular, since one should take into account ultrasoft regions 
(see in Section 1). However, for these cases (when some masses vanish) exact results are 
available in refs. 
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7. Conclusions 



Let us briefly summarize the main results of this paper. 

We have considered the application of prescriptions of ref. Jl[] to the construction of 
the threshold expansion of the sunset diagram (Fig. 1) at the three-particle threshold. 
We have treated the most complicated case, when all particles involved in the threshold 
cut are massive. The expansion involves two contributions: the (h-h) and the (p-p) ones. 

The main technical difficulties are related to the (h-h) contribution. Basically, it 
corresponds to a formal Taylor expansion around the threshold. To construct the coef- 
ficients of this expansion, one needs to calculate threshold values of the corresponding 
sunset integrals (0) with higher powers of the propagators. Although results for the low- 



est integrals were known up to the finite parts in e the recursive procedure (based 



on the techniques p2|,|23|) required either e-part of those integrals or an explicit result 
for one of the higher integrals (e.g. L<r(4 — 2e; 2, 2, 1)). We have calculated this integral 
analytically, in terms of dilogarithms and Clausen functions (j38[) . Note that it contains a 
1/s pole which is associated with the threshold singularity (treated in the framework of 
dimensional regularization) . The result fl3~8|) has completed the algorithm for calculating 
higher terms of the (h-h) contribution. 

Collecting the (h-h) and (p-p) contributions, we get the terms of the threshold expan- 
sion. A comparison with the results of semi-numerical program |1| has been performed. 
It was shown that the combined expansion (using the variables and y/k 2 below 

and above the threshold, respectively) provides good analytical approximations in a 
wide region of values of the external momentum squared. This illustrates how efficiently 
the general procedure of threshold expansion jl| works for diagrams with three-particle 
thresholds in the totally massive case. 
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work has been done, during our visits supported by the EU INTAS grant 93-0744. We 
would like to thank O.V. Tarasov and J.B. Tausk for useful comments. A. D.'s research 
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No. 1/73611. Both of us acknowledge partial support from the Russian Foundation for 
Basic Research, project 98-02-16981. 

Appendix A: Recurrence relations for different masses 

The main parameter of recursion is the sum of the indices, a = v\ + v 2 + ^3- Under 
"lower integrals" we understand the integrals with lower values of a. When one of the 
z/'s is zero, the corresponding integral is trivial (see eq. (0)). The integrals with a = 3 



(ui = v-i = U3 = 1) and a = 4 (yi = v% = 1, 1/3 = 2 and permutations) are known [g4 
We are therefore interested in calculation of the integrals with a > 5. 

Using eqs. fl24D , fl25|) we can express any integral Lt(ti; u x , is 2 , ^3) with positive z/s 



in terms of L T (n; 1, 1, a — 2), L T (n; 2, 1, a — 3), L T (n; 1, 2, a — 3) and L T (n; 2, 2, a — 4). 
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Furthermore, using the integration-by-parts identities we can express Lt{u; 2, 1, a — 3) 
and L T {n; 1, 2, a — 3) in terms of L T {n; 1, 1, a — 2), L T {n; 2,2, a — 4) and lower integrals, 
via 

2(z/ 3 - l)m 3 (mi + m 2 )(2mim 2 + mim 3 + m 2 m 3 )L T (n; 2, 1, z/ 3 ) 
= 2m\(mi + m 2 )(mi + m 3 )L T (n; 2, 2, v 3 - 1) 

+2m^z/ 3 (z/ 3 - l)(mi + m 2 )(m 2 + m 3 )L T (n; 1, 1, z/ 3 + 1) 

-(zaj - l)(2n - 2v 3 - 3) (mi + m 2 )(m 2 + m 3 )L T {n; 1, 1, z/ 3 ) 

-(n - z/ 3 - l)(mim 2 + mim 3 + m 2 m 3 )L T (n; 2, 1, z/ 3 - 1) 

-(n -u 3 - \)m\L T {n; 1, 2, v 3 - 1) + (v 3 - l)m|L(n; 0, 2, z/ 3 ) 

+ (z/ 3 - l)(mim 2 + mim 3 + m 2 m 3 )L(n; 2, 0, z/ 3 ), (60) 

and an analogous equation for 1, 2, z/ 3 ) (only mi and m 2 , together with the argu- 

ments of the integrals corresponding to v\ and v 2 , are to be permuted). 



Then, considering equations ( p2|) (which follow from ref. [23]), we obtain the following 
solution for the remaining integrals (with V\ = v 2 = 1 and V\ — i/ 2 — 2): 

32m 3 (z/ 3 -l)(z/ 3 -2)(n-z/ 3 -l)(mi + m 2 )(m 2 +m 3 )(m 3 +m 1 )/i T L T (ra; 1, 1, u 3 ) 

= 2(is 3 -2)(m 1 + m 2 ) (3ra-4z/ 3 -l)(3n-4z/ 3 + l)(m 1 + m 3 )(m 2 + m 3 )/i T 

+2(n - 3)(n - z/ 3 - l)m 3 [i T ([i T + m 3 ) 

+(n — 3)(n — 2z/ 3 + l)m 3 (mi + m 3 )(m 2 + m 3 ) 



-2(mi + m 2 ) 



L T (n; 1,1,1/3-1) 



'n - 2z/ 3 + l)(mi + m 2 )(m 2 + m 3 )(m 3 + m x ) 



+2(n - z/ 3 - l)mim 2 (/i T + m 3 ) L T (n; 2, 2, z/ 3 - 3) 

+(mi + m 2 ) j(n - 3) (mi + m 3 )(m 2 + m 3 ) [(2n - 11) (mi + m 2 ) - 2m 3 ] 
+(n — v 3 — l)(n + 2z/ 3 — 9)mim 2 m 3 



+ (n - 1/3 - l)/i T 



(7n - 20z/ 3 + 39) (mi + m 3 )(m 2 + m 3 ) 



— m 3 ^2(n - 3z/ 3 + 6) (mi + m 2 ) + (3n - 14z/ 3 + 33)m 3 
x [L T (n; 2, 1, z/ 3 - 2) + L T (n; 1, 2, z/ 3 - 2)] 

+(n - z/ 3 )(mi - m 2 K - (n — 3) (mi + m 2 )(m 1 + m 3 )(m 2 + m 3 ) 
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+2{n — i/ 3 — 1) ^ T {2m 1 m 2 + 3m!m 3 + 3m 2 m 3 + m 3 ) — m 1 m 2 m 3 | 
x [Lr(n; 2, 1, z/ 3 - 2) - L T (n; 1, 2, z/ 3 - 2)] 

—2(i/ 3 - 2) (mi + m 3 )|(mi + m 2 )(m 2 + m 3 ) [rai(n - 5) - m 3 (n - 3)] 
+2m 3 (n - 1/3 - 1)[// T (mi + m 2 ) +mi(m 2 + m 3 )] j>L(n; 2, 0, z/ 3 - 1) 
—2(1/3 - 2)(m 2 + m 3 )|(mi + m 2 )(m 1 + m 3 ) [m 2 (n - 5) - m 3 (n - 3)] 
+2m 3 (n-v 3 -l)[(m 1 + m 2 )/i T + m 2 (m 1 + m 3 )]j>L(n; 0, 2, v 3 - 1), 

32(n — 1/3 — 3)m 1 m 2 (m 1 + m 2 )(m 2 + m 3 )(m 3 + mi)yU T L T (n; 2, 2, z/ 3 ) 
= 2(n - 3)i/ 3 (toi + m 2 )|(mi + m 3 )(m 2 + m 3 ) [(n - 5) (mi + m 2 ) - 2m< 

-2m 3 (n - v 3 — 3) /x T - mim 2 |L r (n; 1, 1, 1/3 + 1) 



+2(mi + m 2 ) 



(ra - 3) (mi + m 2 )(m 2 + m 3 )(m 3 + mi) 

Lr(n; 2, 2, ia, - 1) 



+2(n — 1/3 — 3)/xr(2mim 2 + mim 3 + m 2 m 3 ) 

+ (mi + m 2 ) j(n - 3) (mi + m 3 )(m 2 + m 3 ) [(2n - 11) (mi + m 2 ) - 2m 3 ] 
+ (n — 1/3 — 3)(n + 2z/ 3 - 5)m 1 m 2 m 3 

+(n — v 3 — 3) /x T - 3(3n - 4z/ 3 - 5) (mi + m 3 )(m 2 + m 3 ) 

+m 3 ^2(n - 3i/ 3 )(mi + m 2 ) + 5(rx - 2v 3 - l)m 3 
x [L T (n;2, 1, v 3 ) + L T (n; 1,2, v 3 )} 

+ {n - v 3 - 2) (mi - m 2 ) j - (ra - 3) (mi + m 2 )(m 1 + m 3 )(m 2 + m 3 ) 

+2(rx — i/ 3 — 3) [-mim 2 m 3 + fi T (2m 1 m 2 + m 3 /x T )] | 
x [L T (n; 2, 1, i/ 3 ) - L T (n; 1, 2, i/ 3 )] 

+2z/ 3 (mi + m 3 )|(mi + m 2 )(m 2 + m 3 ) [(n - 3)m 3 - (n - 5)mi] 
+2m 3 (rx - 1/3 - 3)(m^ + m 2 /x T ) }-£("<; 2, 0, 1/3 + 1) 
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+2z/ 3 (m 2 + m 3 )Umi + m 2 )(m 1 + m 3 ) [(n - 3)m 3 — (n — 5)m 2 ] 

+2m 3 (n - i/ 3 - 3)(m£ + m^Jl^; 0, 2, ia, + 1). (62) 

We note that the factors (n — z/ 3 — 1) (on the l.h.s. of eq. (|6~1~D) and (n — z/ 3 — 3) (on 
the l.h.s. of eq. (|B"2"D) yield (n — 4) at z/ 3 = 3 and z/ 3 = 1, respectively. This illustrates 
the problem arising when getting the results for L T (n; 2, 2, 1) and L T (n; 1, 1, 3) (see in 
section 3). 

Appendix B: Threshold results for some other integrals 

Using eqs. fl28|) , (|30"D and (f|l])> we obtain 



Lr(4-2e;l,l,3) 



7T 



4-2e p2 



F (l+e 



{™ 3 - 4£ 


" 1 











mi 

m 3 



m 2 
m 3 



[m 1 m 2 m 3 



,l/2-3e 



24 /2 ' 5e 



7T 



— + 4vr + tt£ + C 



} + (63) 



Here and below, £ = £(mi, m 2 , m 3 ) and C = C(6i, 6* 2 , 6* 3 ), see eqs. (|4T)D and (p9|). 

Employing explicit recurrence relations from Appendix A, we can obtain results for 
the integrals Lt with higher values of a = v\ + z/ 2 + z/ 3 . For instance, for o = 6 we get 



L / ( I — 2;: 2.2.2) = - 7r 4 ~ 2£ T 2 (l + e) {m^m^- 1 ' 36 /i T 1+5e 



(mi + m 2 )(m 2 + m 3 )(m 3 + mi) + 4mim 2 m 3 



4y / 'mim 2 m 3 /i| 1 



7T 

2e 



3tt + ?r£ + C 



/mim 2 m 3 



/4 



3 + 



mim 2 m 3 + //^ 



(mi + m 2 )(m 2 + m 3 )(m 3 + mi) 



+ 1 



m 3 (m 1 - m 2 ) ^ mi mi(m 2 — m 3 ) 



ln^ + 



m 2 (m 3 - mi) m 3 
in ■ 



(m 2 +m 3 )(m 3 +mi) m 2 (m 3 +mi)(mi+m 2 ) m 3 (mi+m 2 )(m 2 +m 3 ) mi 
1 



2^ T 

L T (4-2e;l,2,3) 



(mi — m 2 ) In — - + (m 2 — m 3 ) In — - + (m 3 — mi) In — - 

m 2 m 3 mi 



(64) 



1 



3e, -l+5e 



-n A - 2e T 2 (l + e) (mim 2 m 3 )-^ 
5m 2 m 3 



x ■ 



3(mi + m 2 )(m 2 + m 3 )(m 3 + mi) - 4mim 3 /i T 



4Wmim 2 m 3 /ij 1 



— + - + tt£ + C 

2 £ 3 



7r mi 



1 /m 1 m 2 m 3 /i T [3 



2. . (mi+m 2 )(mi+m 3 )-m 2 +m 3 

-(3m 2 -m 3 ) - m 2 m 3 



(mi + m 2 )(m 2 + m 3 )(m 3 +mi) 
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mi(m 2 — m 3 ) 
(mi + m 2 )(m- i + m 1 

3 



ln^ + 

m 2 



7713(7712 — mi) 



In 



m 3 



/ \ 1 m l / m 3 

(m 2 + 7713) m h (mi + m 2 ) In — 

ni2 ni2 



(m 2 + m 3 )(m 3 + mi) m 2 
+ 0(e). (65) 



The result for Ly(4 — 2e; 2, 2, 2), eq. (|64"1) , is totally symmetric in mi, m 2 , m 3 , as it should. 
Using eq. (|65| ) (and its permutations) and the results for lower integrals, one can easily 
obtain result for L^(4 — 2e; 1, 1,4) (and its permutations) via eqs. fl24D -fl25D. 
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gure 1 : The sunset diagram 




gure 2: Threshold approximations for unequal masses (£1 = 0.1, £2 = 0.3, £ 




Figure 3: Threshold approximations in the case of equal masses 



